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Abstract: We study tunneling of massless and massive particles through the smeared 
quantum horizon of the extra-dimensional Schwarzschild black holes. The emission rate 
of the particles' tunneling is modified by noncommutativity effects in a bulk spacetime 
of dimension d. The issues of information loss and possible correlations between emitted 
particles are discussed. We show that even by considering both noncommutativity and 
braneworld effects, there is no correlation between different modes of evaporation at least 
at late-time and within approximations used in the calculations. However, incorporation 
of quantum gravity effects such as modification of the standard dispersion relation or 
generalization of the Heisenberg uncertainty principle, leads to the correlation between 
emitted particles. Although time-evolution of these correlations is not trivial, a part of 
information coming out of the black hole can be preserved in these correlations. On the 
other hand, as a well-known result of spacetime noncommutativity, a part of information 
may be preserved in a stable black hole remnant. 
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1. Introduction 

Radiation spectrum of an evaporating black hole which is almost similar to black body 
radiation spectrum, can be described by a characteristic temperature dubbed Hawking 
temperature. By adopting a quantum field theoretical approach in curved spacetime, this 
characteristic temperature is given by Th = <^bG ' w here k is the surface gravity that 
demonstrates the strength of the gravitational field near the black hole surface [0 . Unfor- 
tunately, Hawking's approach leads to a non-unitarity of quantum theory which maps a 
pure state to a mixed state as a result of purely thermal nature of the spectrum. A few 
years ago, Parikh and Wilczek ||] suggested a new method based on null-geodesics to ex- 
tract the Hawking radiation via tunneling across the event horizon. This method modifies 
the form of the black hole radiation spectrum that is due to inclusion of back-reaction ef- 
fects. The tunneling process clarifies that the modified radiation spectrum is not precisely 
thermal and this leads to the unitarity of underlying quantum theory ||]. Nevertheless, the 
form of the correction is not sufficient by itself to recover information because of its failure 
to have correlations between the tunneling probability of different modes in the black hole 
radiation spectrum Recently, Nicolini, Smailagic and Spallucci (NSS) ||, in a new 
conceptual approach to space noncommutativity, have shown that black hole evaporation 
process should be terminated when black hole reaches to a minimal mass, i.e. black hole 
remnant (for a review see ||). This minimal mass is a consequence of existence of minimal 
observable length P, || : ^, [10|, O]. The NSS viewpoint on coordinate noncommutativity is 
carried out by the Gaussian distribution of coherent states and is consistent with Lorentz 
invariance, unitarity and UV-f initeness of quantum field theory fl^l . Moreover, non- 
commutativity of spacetime is an inherent property of the manifold by itself even if there 
is no gravity. In this framework, some kind of divergencies which arise in general relativity 
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and black hole physics, can be removed also. In a recent paper [13], we have studied the 
Parikh-Wilczek tunneling through the quantum horizon of a Schwarzschild black hole in 
noncommutative (NC) spacetime. We have shown that, in the proposed NC setup there is 
no correlation between different modes of the radiation which reflects the fact that informa- 
tion does not come out continuously during the evaporation process, at least at late-time. 
The fundamental motivation of the present paper is to find an alternative framework to 
recover the lost information in the black hole evaporation process (see []14| for comprehen- 
sive reviews on existing approaches to resolution of black hole information paradox). This 
alternative framework can be identified with the correlations between different modes of 
evaporation. Another question that we focus on is the possible impact of braneworld sce- 
narios and correlation between different emitted modes of black hole evaporation. In other 
words, we want to know whether inclusion of extra spatial dimensions in the framework of 
extra dimensional scenarios can create correlation between different modes or not. This is 
an important question which is considered here via the tunneling framework of Parikh and 
Wilczek. With these motivations, we apply back-reaction and noncommutativity effects to 
study black hole evaporation process in a bulk spacetime of dimension d. Recently, we have 
shown that incorporating quantum gravity effects (such as Generalized Uncertainty Prin- 
ciple (GUP) @, |, H |H| |ll| (see also ||, H) and Modified Dispersion Relations (MDR) 
fl7| , |i~8| , [l9|| ) in the black hole evaporation process has the capability to produce correlations 
between different emitted modes [2C]. Although the time evolution of these correlations 
are not obvious and deserve further investigations, the possibility of realizing these corre- 
lations itself is an important step toward the resolution of black hole information paradox. 
Here we generalize this achievement to our extra dimensional framework and discuss some 
related issues. 

The paper is organized as follows: In Section || we construct a noncommutative frame- 
work in the presence of extra dimensions and we study evaporation of Schwarzschild black 
hole in this setup. In Section ||, a detailed calculation of quantum tunneling near the 
smeared quantum horizon is provided. We calculate tunneling probability of massless and 
massive particles through the quantum horizon. The problem of lost information and possi- 
ble correlations between emitted particles are investigated. The capability of the black hole 
evaporation process to produce the correlations between different emitted modes via incor- 
poration of quantum gravity effects is investigated in Section Summary and conclusions 
are presented in Section |B[ 



2. Noncommutativity in the Presence of Extra Dimensions 

Noncommutative quantum field theories (NCQFTs) have attracted a lot of attentions dur- 



ing last two decades pl[ (see also [22] for a purely phenomenological viewpoint on this 



issue). There are also a lot of attempts to possible realization of these effects on experimen- 



tal grounds [23]. The NC spacetime coordinates, in their simplest form, are deformations 
of usual spacetime coordinates (with an arbitrary number, d of spacetime dimensions) in 
which operators X A representable by the Hermitian operators X A , characterize a nonvan- 



- 2 - 



ishing commutation relation as follows 

riAB 

[X A ,X B ]=iO AB = i^, (2.1) 

A NC 

where A^vc is a characteristic energy or inverse length scale related to the NC effects 
and C AB (6 AB ) is a dimensionless, real and anti-symmetric matrix, whose elements are 
constant and are assumed to be of the order of one. The physical explanation of the 9 AB is 
the smallest fundamental cell of observable area in the ^4I?-plane, in the same way as Planck 
constant h explains the smallest fundamental cell of observable phase space in quantum 
mechanics. The scale A^c, possibly and most reasonably is of the order of Planck scale 
Mp\. This is supported by the fact that most of the phenomenological studies of the NC 
models have presumed that A^rc cannot lie far above the TeV scale |2^, 23 1. Since the 
fundamental Planck scale in models with large extra dimensions (LEDs) becomes as small 
as a TeV in order to solve the hierarchy problem [24, |2[|; it is possible to set the NC effects 
in a TeV regime. 

The framework of NCQFT is achieved by substitution of ordinary products in com- 
mutative theory with Moyal ^-products [26|. There exist many formulations of NCQFT 
based on the Moyal ^-product that lead to failure in resolving of some important problems, 
such as Lorentz invariance breaking, loss of unitarity and UV divergences of quantum field 
theory. Unfortunately, no flawless and completely convincing theory of noncommutativity 
exists yet. Recently, the authors in Ref. [^] explained a fascinating model of noncommuta- 
tivity, the coordinate coherent states approach, that can be free from the problems 
mentioned above. In this approach, General Relativity in its usual commutative form as 
described by the Einstein-Hilbert action remains applicable. If noncommutativity effects 
can be treated in a perturbative manner, then this is defensible, at least to a good approx- 
imation. Indeed, the authors in Ref. 27] have shown that the leading noncommutativity 
corrections to the form of the Einstein-Hilbert action are at least second order in the non- 
commutativity parameter 9. The generalization of the QFT by noncommutativity based on 
coordinate coherent state formalism is also interestingly curing the short distance behavior 
of pointlike structures || [l^] (see also [28, |2{|). In this method, the particle mass M, in- 
stead of being quite localized at a point, is described by a smeared structure throughout a 
region of linear size \J~Q. In other words, we shall smear the point mass Dirac-delta function 
utilizing a Gaussian function of finite width. So, instead of a delta function distribution, 
the mass density distribution of a static, spherically symmetric, particle-like gravitational 
source in ci-dimensional bulk spacetime will be given by a Gaussian distribution of minimal 
width Vd as follows 1 @ 

M r 2 

Pe(r) = - -rre-^. (2.2) 

(4vr6») — 

The line element which solves Einstein's field equations Gba = ^G^Tba i n the presence 
of smeared mass sources can be obtained as 

* 2 — (l - ^fV 2 + (l - r*' 2 + r^ d _ 2) , (2.3) 



throughout the rest of this paper, natural units are used so that H = c = fcs = 1. 
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with 

M e = —r-r- / dt eT l t^—> , m = — £ . (2.4) 

Where dQ^ d _ 2 ^ is the line element on the (d — 2)-dimensional unit sphere and d is spacetime 
dimensionality. The volume of the (d — 2)-dimensional unit sphere is given by 

d-l 



Gd is gravitational constant in d-dimensional spacetime which in ADD 25] model is given 

by 

G d = ^^M*jd, (2.6) 

where Mpi is the (i-dimensional Planck mass and there is an effective 4-dimensional Newton 
constant related to Mp\ by 

MpJ d = 4vrG 4 i? d - 4 , (2.7) 

where R is the common radius of all (d — 4)-toroidally compactified extra dimensions. We 
note that in this work, the conventions for definition of the fundamental Planck scale Mpi 
are the same as which have been used by ADD (and also by authors of Ref. pOfl ). The 
metric (^^) smoothly interpolates between de Sitter core around the origin and an ordinary 
Schwarzschild geometry far away from the origin 0, |29| . On the other hand, the curvature 
singularity at the point r = is eliminated by noncommutativity as an inherent property of 
the manifold. In this situation, a regular de Sitter vacuum state will be formed accounting 
for the effect of the NC coordinates fluctuations at short distances. Also, this scenario tends 
to usual Schwarzschild spacetime at large distances. In the limit 9^0, one recovers the 

complete Gamma function, f 46 dte~ t t^^~^ — ► T(^p), and M.g — ► m as expected. Taking 
a 4-dimensional viewpoint in the commutative case, one can simply obtain Mq = G4M. 
Then the NC (modified) Schwarzschild solution reduces to the commutative (ordinary) case. 
Depending on the different values of M, d and Mpi and within a numerical procedure, the 
metric (2.3) displays three possible causal structure |, 0, ||, (see also ||, ||]): 1- 



It is possible to have two distinct horizons (Non-extremal Black Hole), 2- It is possible to 
have one degenerate horizon (Extremal Black Hole), and finally 3- It is impossible to have 
horizon at all. 

The radiating behavior of such a modified Schwarzschild black hole can now be easily 
investigated by plotting the temporal component of the metric, goo , versus the radius r 
for an extremal black hole with different values of M and d. This has been shown in 
Fig. (||). As this figure shows, the coordinate noncommutativity leads to the existence 
of a minimal non-zero mass in which black hole (due to Hawking radiation) can shrink 
to 2 . It should be emphasized that based on our computations, if we set d = 11, we find 



2 For simplicity of numerical calculations, we assume 6 = 1. Naturally (Ajvo) points out the energy 
threshold beyond which a particle moves into the deformed spacetime in which quantum gravity effects 
become important. 



- 4 - 




Figure 1: goo versus the radius r for different values of black hole mass M and different number of spacetime 
dimensions, d. The figure shows the possibility of having extremal configuration with one degenerate event horizon 
(Extremal Black Hole) at M = Mq. This shows the existence of a minimal non-zero mass that black hole shrinks to. 
On the right-hand side of the figure, curves are marked from bottom to top by d = 4 to d = 10 (with M Pl = 1). 



Table 1: Minimal non-zero mass of the black hole (remnant mass) for different number of spacetime 
dimensions. 



Extremal Black Hole ; M Pi = 1 TeV 


d 


Minimal Mass (TeV) 


4 


M = M S3 1.9 


5 


M = M ~ 3.94 


6 


M = Mo S3 5.43 


7 


M = Mo sa 5.77 


8 


M = Mq S3 5.08 


9 


M = Mq S3 3.88 


10 


M = Mq S3 2.65 


11 


M = Mq S3 1.64 



M = Mo ~ 1.64 TeV, i.e. the minimum value of black hole mass reduces to a value less 
than its value for d = 4 (see Table |]). Therefore, in theories with large extra dimensions, 
if the number of spacetime dimensions becomes sufficiently great, e.g. d > 11 with a 
sufficiently small NC parameter (Ajvc* ~ Mpi ~ 1 TeV), then the noncommutativity effect 
can enhance the possible formation and detection of black holes in TeV-scale collisions at 



the LHC 3 p4] 



3 A micro-black hole in theories with large extra dimensions can be produced at the large hadron collider 
(LHC) just in the situation that E cm > Mo, where E cm is parton-parton center-of-mass energy which is 
equal to 14 TeV. 
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The event horizon radius, r#, can be obtained from the equation goo (r# ) = 0, which 



gives 



with 



-i 2.M e (r H ) _ 



.2 \ 4 



' II . 



d+l d + 3 rjj 
4 ' 4 '4? 



d- 1 



2 (l + d) (d + 3) 



W 



d + 5 d + 3 
4 ' 4 '40 



+ 



d- 1 



26 



+ l + d 



d+l 

r 2 H \ 4 rf+i _r| 



(1 + d) 



-l 



(d-2)vr— 



Where W(/i, i^, z) shows the Whittaker function defined in terms of the hypergeometric 
functions as follows 



1 



v, z) = hyper geomi v — /i + - , 2v + 1 , z) e 2 z u+ " 



Analytical solution of the Eq. (2.8) for r# in a closed form is impossible, but we can 
approximately solve it by setting r# = 2M into the lower incomplete Gamma function of 
the Eq. (EO) to find 4 



r H = (2M e (r H ) 



ih \87rG d M a 



d-3 
) 2 



M — 2- _m+ a ,/d+l d + 3 M 



d+ 
4 



— e 2e yy 



d- 1 



;l +<1) (<! + :]) ) + ( 2^- + 1 + d 



M 



/-2 



e 29 



(d-2)n—6 



-1 



1 

d-3 



(2.c 



For very large masses, one can easily recover the classical Schwarzschild radius in a space- 

1 

time with d dimensions, r# ~ (2m) d ~ 3 . 

Now we calculate Hawking temperature of noncommutative black hole by investigating 
its radiating behavior. By definition, 



1 dgoo, 
47r dr 



\r=ru ■ 



(2.10) 



Black hole temperature for some arbitrary number of spacetime dimensions can be calcu- 
lated as follows 



T 4 = G A M 



£ \2VeJ r H e~ 



2irr 



H 



29 



(2.11) 



4 We will use this approximation in the next section too. 
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T 5 = G 5 M 



r H e 



LtL 

4,0 



3ir 2 i 



H 



n = GqM 



or is. 

_ U/g 

8tt% 



I2(n6 



r H e 



40 



1 + -s- + 

4 



J3 



r 2 
49 



32(vr6^ 



1 + 



69 369' 



+ 



H 



G 7 M 



32 



f 2 

LH 



IOtt 3 ^ 



1 + 



649 2 2566> 3 \ 



H 



H 



8O(vr0) 3 

and so on. £{x) shows the Gauss Error Function defined as follows 

s/k Jo 



H 



(2.12) 



(2.13) 



(2.14) 



- t2 dt. 



While for all even values of d Hawking temperature can be expressed in terms of the Gauss 
error function, for odd values of d Hawking temperature has a closed analytical form. For 
very large distances (or commutative case), the function £(tt^) fo r even values of d tends 
to unity and other terms on the right hand side will be vanishing exponentially. Therefore, 
one recovers the classical Hawking temperature in d dimensions 



d-3 
d-2 



TT 



d- 1 



2-d 
n H ■ 



(2.15) 



In this situation, if we want to obtain the simple (i-dependent form of the noncommutative 
Hawking temperature, then it can be approximated by utilizing the Eq. (2.9) as follows 

d-3 



4nrH 



(2.16) 



We will come back to this relation later in this paper. 

It should be stressed that, if instead of distribution Q2.2] ), we adopt other kinds of 
probability distributions, only the smeared mass distribution Mq will change and all general 
properties will be the same as above. For example, we consider a Lorentzian distribution 
of smeared matter as follows 5 



MVW 



P9'(r) = 

The smeared mass distribution is now given by 



(2.17) 



Mf)> 



2m(d-2)!! [VF t 



j.d-2 



7r(d-3)!! J (l + t 2 ) 



dt. 



(2.18) 



In the limit of 9' — ► 0, we obtain Mg> — > m. The Lorentzian smeared mass, has 
limiting properties similar to Gaussian one and therefore can be compared with Gaussian 
smeared mass distribution, M.g. Thus, most of the outcomes that we obtained for Gaussian 



5 Generally, the noncommutativity parameter 9' appeared in Eq. (2.17) is not exactly the same as 9 
appeared in Eq. (2.2). 
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profile, stay applicable if we pick out a different form of the probability distribution of 
smeared matter 6 (see also [13, |29[|). 

As an important remark, we note that some authors have studied black hole thermo- 
dynamics in the NC framework adopting a coordinate noncommutativity against coherent 
state approach (see |35| ] and references therein). A question then arises: what is the dif- 
ference between these two approaches? The standard way to handle NC problems is based 
on Moyal ^-product. That means to use complex number commuting coordinates and shift 
noncommutativity in the product between functions. This is mathematically valid, but it 
is physically useless since any model written in terms of ^-product, even the simplest field 
theory, is nonlocal and it is not obvious how to handle nonlocal quantum field theory. One 
suggested approach is perturbation in the 9 parameter |3(|. This is physically reasonable 
due to the fact that once expanded up to a given order in 9, the resulting field theory be- 
comes local. The smeared picture of particles based on coordinate coherent states defines 
complex number coordinates as quantum mean values of the original noncommuting ones 
between coordinate coherent states |l2|]. In other words, in this setup one can see commut- 
ing coordinates as classical limit (in the quantum mechanical sense) of the noncommuting 
ones. In this framework, the emergent semiclassical geometry keeps memory of its ori- 
gin. For example, free propagation of a point-like object is described by a minimal width 
Gaussian wave-packet as has been considered in our setup. So, the difference between two 
approaches lies in the definition of quantum field theoretical propagators. 



3. Particles' Tunneling Near the Smeared Quantum Horizon 

Now we provide a detailed calculation of quantum tunneling from the NC black hole event 
horizon. To portray the NC quantum tunneling procedure, where a particle moves in 
dynamical geometry and passes through the horizon without singularity on the path, we 
should take advantage of a coordinates system that, unlike Schwarzschild coordinates, is 
not singular at the horizon. A convenient choice in this regard is Painleve coordinate 
transformation [[37]] which is obtained by definition of a new NC time coordinate 



dt = dt * + r d-3 _ 2Mg dr > 

where t s is the Schwarzschild time coordinate. Note that only Schwarzschild time coordi- 
nate is transformed, while radial coordinate and angular coordinates remain unchanged. 
Now the NC Painleve-Schwarzschild metric in spacetime with d dimensions takes the fol- 
lowing form 



-(i-^V+V^^ + * a+ ^«>- (3 - 2) 



6 The difference between various forms of smeared mass distributions is considerable around M < Mo 
where there are sensitivity to noncommutative effects and detailed form of the matter distribution. There- 
fore, we focus on situations that the condition M > Mo holds. 
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The metric in these new coordinates is stationary and non-static. In addition, neither 
coordinate nor intrinsic singularities are present. Radial null geodesies (the equation of 
motion for massless particles) in this geometry are as follows 



where the upper (lower) sign corresponds to an outgoing (ingoing) null geodesic respec- 
tively. If we assume t increases in the direction of future, then the above equations will 
be corrected by the particle's self-gravitation effect. Kraus and Wilczek |58[ scrutinized 
the motion of particles in the s-wave as spherical massless shells in dynamical geometry 
and expanded self-gravitating shells in Hamiltonian gravity (see also |39t| ). Supplementary 
elaborations was accomplished by Parikh and Wilczek In this manner, we are going 
to develop their framework to the NC coordinate coherent states in a model universe with 
extra-dimensions . 

We set the total ADM mass (M) of the spacetime to be fixed, and permit the hole 
mass to fluctuate. Also, we consider the reaction of the background geometry to an emitted 
quantum of energy E which moves on the geodesies of a spacetime with M replaced by 



M — E. In this manner, we should replace M by M — E in Eqs. (3J, 3^, and |3.3j ). The 
geodesic equation of motion, Eq. (|3,3|) , can be written as 



^ ±1 -^t^> (3 ,) 

Since the characteristic wavelength of the radiation near the horizon is always small due to 
infinite blue-shift, the wave-number approaches infinity. Hence the WKB approximation 
is valid near the horizon. In the WKB limit, the tunneling probability or emission rate for 
the classically forbidden trajectories as a function of imaginary part of the particle action 
at stationary phase takes the following form 

T ~ exp(-2Im I). (3.5) 

To compute the imaginary part of the action, we consider a spherical shell of massless par- 
ticles that move on radial null geodesies. Since the radial null geodesies can be interpreted 
as s-waves outgoing positive energy particles which cross the horizon outward from rj n to 
r ut, we use these radial null geodesies to calculate the Im/, as follows 7 

rr ou t rr ut rPr 

Im/ = Im / p r dr = Im / / dp' r dr. (3.6) 

J Ti n J Ti n JO 

Therefore, we incorporate back-reaction effects in a finite separation between initial and 
final radius as a result of self-gravitation effects of outgoing shells. On the other hand, 
according to energy conservation, the tunneling barrier is produced by a change in the 



7 Note that we require ri„ > r ut where n n = (d-tyn*™ 2) r(^i) Jo * * ( 2 > and r out = 

(d- 2 )n (d _ 2) TiSET) Jo 
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radius (reduction of black hole horizon) just by emitting particles. We now change the 
integration variable from momentum to energy by applying Hamilton's equation of motion, 
* = & • So, we find 

pM-E pTout A T 

ImI = Im / / — dH, (3.7) 

JM Jr in r 

where the Hamiltonian is given by H = M — E' . The r-integral can be solved firstly by 
deforming the contour so that 

Im/ = Im / / (-dE?). (3.8) 

Jo Jr in i __ . hMe(M-E') 



The r-integral has a pole at the horizon in which it lies on the line of integration and leads 
to (— m) times the residue 

ImJ = Im (^3) (zM e (M-E'))~\-dE'). (3.9) 

The above integration can be preformed without writing out the explicit form of radial null 
geodesies. In the vicinity of the horizon, r behaves as 

f~(r-r H )K(M)+0({r-r H ) 2 ), (3.10) 

where n{M) is the horizon surface gravity. The mentioned integration can be done simply 
by inserting the relation ( 3.10| ), in which self-gravitation is included, into the integral ( |3.7| ) 
and performing the r-integral as follows 

f E f r °* drdE' f E dE' 
Im/ = -Im/ / — — =7t/ — — — . (3.11) 

Jo Jr in (r - r H ) K (M - E>) J «(M - E>) V ; 

Comparing two integrands in the relations (^) and ( |3.11| ) gives us 

k(M-E) = d —^- -. (3.12) 

2{2M e {M - E))~ 

Now, Hawking temperature is given by 

Th = am „ — i^^. ( 3.i3) 

27r 4ir(2M e (M))< 



I d-3 



It is clear that, two Eqs. ( p. 16 ) and ( |3.13j ) for the NC black hole temperature coincide, 
but note that the approach to achieve ( |3.13j ) uses some approximations. This coincidence 
confirms that our calculations are approximately correct. 

Here, we would like to obtain the radial geodesies of the massive particles which are 
different with massless case. Recently, some authors have extended Parikh-Wilczek work 
to the massive particles' tunneling |4(], [II], fi~2^ . In this regard, once again we follow the 
line element of the NC Painleve-Schwarzschild black hole, that is 

ds 2 = g 00 dt 2 + 2g 01 dtdr + g n dr 2 + r 2 dtt 2 d _ 2) . (3.14) 
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This line-element can be obtained by definition of a noncommutative time coordinate trans- 
formation as 

dt = dt s + dt syn , (3.15) 

where 



y/2M g r d - 



3 



5oi 



dtsyn = r^-2M e dr = 'm^ ^ 
This relation, which has been written based on the Landau's theory of the synchronization 
of clocks p3] , allows us to synchronize clocks in any infinitesimal radial positions of particles 
(dfJ((2_2) = 0). Since the tunneling phenomena through the quantum horizon, the barrier, 
is an instantaneous procedure, it is important to apply Landau's theory of the coordinate 
clock synchronization to the tunneling process. The mechanism for tunneling through the 
quantum horizon is that particle anti-particle pair is created at the event horizon. So, 
we have two events that occur simultaneously; one event is anti-particle and tunnels into 
the barrier but the other particle tunnels out the barrier. In fact, the relation ( 3.1€| ) 



mentions the difference of coordinate times for these two simultaneous events occurring at 
infinitesimally adjacent radial positions. According to the non-relativistic quantum theory, 
de Broglie's hypothesis and the WKB approximation, it can be easily shown that the 
treatment of the massive particle's tunneling as a massive shell is approximately derived 
by the phase velocity v p of the de Broglie s-wave. The relationship between phase velocity 
v p and group velocity v g is given by [40, 41, 42] 



v p = r = -v g . (3.17) 



In the case of dQu_2) = 0> according to the relation ( 3.16| ), the group velocity is 

500 

v g = . 

5oi 

Thus, the outgoing motion of the massive particles take the following form 

900 r d ~ 3 - 2Mg 



2501 2^2M e r 



d-3 



(3.18) 



(3.19) 



We should incorporate back-reaction effects in this relation and then the result should be 
inserted into relation ( |3.7| ) to obtain a relation corresponding to ( |3.9| ) but now for massive 
particles. Let us come back to the integral ( |3.9| ). This integral has not a closed analytical 
solution, so we are going to compute Taylor series expansion of integrand, with respect to 
the particle's energy E, about the zero point and just to first order. Then we can solve the 
integral to find 



Iml = 7r 



d-l 



M \ A/ 2 

l6irG d MTr-^[ — ) 2 e-—[AM*H+2deM' l +§eM' z +Ad9 z +'Ze z +d z z )(d-2) 



i 

e- 2 (d-l)~ 1 (H-d)~ 1 (d+3)- 1 d ~ 3 E(-Ed A e 3 -180M 3 6 2 -90M6 3 -4:EM 4 H6d 2 -16EM 6 H'+ 
40 EM 6 7i+2 M6 3 d i +4 M 3 9 2 d 3 -120 M 5 6 H+20 M 3 6 2 d 2 -36 M 3 6 2 d+12 M8 3 d 3 -8 M0 3 d 2 + 



- 11 - 



16M 5 9dH + 8M 5 9 d 2 H + 8 EM % Ud - 80 EM A H9 - 16 Ed0 2 M 2 - 2 Ed 2 6 2 M 2 + 32 £0 dM 4 + 
4 EO d 2 M 4 -108 M9 3 d+60 E6 M 4 -30 E9 2 M 2 -\h EdB 3 -23 Ed 2 9 3 -9 Ed 3 9 3 -36 EM 4 H9 d 

(d-3y 2 9- 1 {5 + dy 1 M- 1 (<iM 4 H + 2d9M 2 + 69M 2 + 4d9 2 + 39 2 + d 2 9 2 ^j , (3.20) 
where H = hypergeom (1, ^±^ ) and W = hypergeom (2, ^f-J. Since, both par- 



ticle and anti-particle (corresponding to a time reversed state that can be seen backward 

in time by replacing y^Pf by — v/ ^^f in the metric ( |3.2[) ) anticipate in the emission 
rate of the Hawking process via tunneling with the same amounts, we have to add their 
amplitudes firstly and then squaring it to obtain the emission probability (Eq. ( |3.5| )). 

On the other side, using the first law of black hole thermodynamics, dM = ^dS, 
and Eq. ( |3.11| ) , one can find the imaginary part of the action as follows [31 |4| (see also 



@, m m, 

-, r S NC (M-E) i 

ImJ = --/ dS = --AS NC , (3.21) 

S*jvc(M) ^ 



2 



which shows that the tunneling approach concurs with the thermal Hawking spectrum 
just in the limit of low energy. Hawking radiation of black holes as quantum tunneling 



was also studied in the context of string theory 44], and it was demonstrated that the 
tunneling probability in the high energy limit is analogous to a difference between counting 
of states in the microcanonical and canonical ensembles. It can be easily seen that our 
computed tunneling probability to first order in E generates the Boltzmann factor in the 
canonical ensemble V ~ exp(— (3E), which is characterized by the inverse temperature as 
the coefficient (3. So, the tunneling probabilities in the high energy limit are proportional 
to exp(A5') (see also @), 

f inal 

T ~ exp(-2Im/) ~ = exp(A,S) = exp[S(M — E) — S(M)], (3.22) 

where AS is the difference in black hole entropies before and after emission. This means 
that at higher energies the emission rate depends on the final and initial number of ac- 
cessible microstates of the system. Therefore at higher energies the emission spectrum 
cannot be exactly thermal since the high energy modifications flow from the physics of 
energy conservation with noncommutativity corrections. In fact, to linear order in E, two 
expressions for T in the microcanonical and canonical ensembles coincide. Thus, obviously 
the emission rate ( 0.22 ) deviates from the pure thermal emission but is corresponding to 



an underlying unitary quantum theory [||, f| . 

At this stage, one can easily show that there are no correlations between emitted 
particles even with incorporation of the noncommutativity modifications in an extra di- 
mensional setup. If there is any correlation between different emitted modes, at least part 
of the information coming out of the black hole to be retained in these correlations. As our 
analysis shows, at least up to second order in particle's energy, there will be no correlation 
between the different modes of radiation. This reflects the fact that information does not 
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come out continuously during the evaporation process at least at late-times 8 . This indi- 
cates that the tunneling probability of two particles of energy E\ and E 2 is comparable to 
the tunneling probability of one particle with composed energy, E = E% + E2, i-e. 

AS El + AS E2 = &S {El+E2) => X (Ei + E 2 ;E 1 ,E 2 ) = 0. (3.23) 

It can be numerically confirmed that these probabilities of emission are actually uncorre- 
cted. On the other hand, the statistical correlation function x(E; E\, E 2 ), is zero which 
leads to the independence between different modes of radiation during the evaporation. Ac- 
cordingly, space noncommutativity and higher dimensional considerations are not enough 
to provide the basis for existence of the correlations between different modes. 

There are four main proposals about what happens to the information that falls into 
a black hole. Based on the first proposal, the black hole can evaporate completely as un- 
corrected thermal radiation in each mode and all the information including the original 
quantum state that formed the black hole (excluding its mass, charge and angular momen- 
tum), would disappear from our universe. However, this proposal allows the pure states to 
evolve into the mixed states, which is incompatible with the basic principles of quantum 
mechanics. Second proposal is that the black hole can completely disappear, but the in- 
formation appears in the final burst of radiation when the black hole shrinks to the Planck 
size. A third possibility is that the information comes out in non-thermal correlations 
within the Hawking radiation, the process being portrayed by a unitary S-matrix. In other 
words, there are non-thermal correlations between different modes of radiation throughout 
the evaporation process that information emerges ceaselessly encoded through them. And 
fourth idea is that the black hole never disappears completely, and the information is not 
lost, but would be stored in a Planckian-size stable remnant. In this section, we have 
studied the credibility of the third and fourth conjectures within a NC spacetime feature. 
We have shown that inclusion of noncommutativity and braneworld effects are not enough 
to create correlations between emitted modes of black hole evaporation. However, due to 
spacetime noncommutativity, the information might be preserved by a stable black hole 
remnant. In our opinion, this is one of the successful scenarios for information loss paradox 
currently [13]. In next section we introduce a procedure to create correlations between 



emitted modes. 



4. Quantum Gravity as a Source of Correlations 

All quantum gravity scenarios support the idea that there is a fundamental length scale that 
cannot be probed. For instance, in string theory there is a restraint in probing distances 
smaller than the string length. Hence the standard uncertainty principle is modified to 
include this restricted resolution of the spacetime points. The result of this modification is 
known as the Generalized Uncertainty Principle (GUP) which is actually a representation 
of the quantum nature of spacetime at Planck scale. In a model universe with extra 



Our calculations are performed based on Eq. (3.5) in stationary phase, hence the phrase late-time. 
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dimensions, GUP can be generalized to the tensorial form as follows 

Axi>— + aL 2 Pl A Ph (4.1) 
Api 

i 

where the fundamental Planck length is defined as Lpi = G^f 2 and a is a dimensionless 
constant of the order of unity that depends on the specific model. The above GUP can be 
acquired in the context of string theory |?J , NC quantum theory || , loop quantum gravity 
0, or from black hole gedanken experiments (TO]. In the standard limit, Axj S> L P i, 
it yields the ordinary uncertainty principle, AxiApi > 1. The second term in r.h.s of 
Eq. becomes considerable when the momentum and distance scales approach the 

Planck scale. For a spherical black hole the thermodynamic quantities can be obtained 
in a heuristic manner using the standard uncertainty principle [49|. Application of the 



GUP to black hole thermodynamics in the same manner, modifies the results dramatically 
by incorporation of quantum gravity effects in the final stages of black hole evaporation 
with a rich phenomenology 16]. In this section we use GUP to find an alternative to 



retrieve the lost information in the black hole evaporation process. This alternative is 
related to the correlations between different modes of evaporation. Can GUP provide 
the required correlations? This is an important question which we answer in a tunneling 



framework. A first effort in this direction was given by Arzano et al [45], resulting to the 
quantum corrected entropy with an additional logarithmic correction term. However, such 
a formulation does not prepare a way to find a non-zero statistical correlation function. 
In this section we follow a different way to retrieve the information. In a recent paper 
|2C], we investigated the modifications of the Hawking radiation by the tunneling process 



and the GUP. By using the GUP-corrected de Broglie wave length, the squeezing of the 
fundamental momentum-space cell, and consequently a GUP-corrected energy, we find 
the non-thermal effects which lead to a non-zero statistical correlation function between 
probabilities of tunneling of two particles with different energies. In this section, we extend 
our framework to model universe with extra dimensions. 

Amelino-Camelia et al [^] (see also studied the black hole evaporation process 
after an analysis of the GUP-induced modification of the black body radiation spectrum. 
If GUP is fundamental concept of quantum gravity, it should appear in de Broglie relation 
as follows 

A ~ - (1 + aL 2 Pl p 2 ) , (4.2) 

or 

£ ~ E(l + aL Pl E 2 ). (4.3) 

There are other compelling reasons from NC geometry and loop quantum gravity that 
support relation (4.3) (see for instance [5C, 51] and references therein). Now, we consider a 



massless particle i.e. a shell and take into consideration the response of the background ge- 
ometry to a radiated quantum of energy E with GUP correction i.e. £. The particle moves 
on the geodesies of a spacetime with M — £ substituted for M. With these preliminaries 
and also considering the deformed Hamilton's equation of motion as 



(l + *L Pl £ 2 )^L\ r , (4.4) 
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eventually the imaginary part of the action takes the following form 



Iml = Im / (—iri, , 

Jo V d - 3 



l + aL 2 Pl £' 2 ) (2(M-£" V ; 



(-d€') 



(4.5) 



The leading order correction is proportional to the square of (y/ot L p ). This integral, similar 
to the integral ( |3.9[ ), has not a closed analytical solution, so we are going to compute Taylor 
series expansion of integrand (with respect to the particle's energy E) about the zero point 
and just to first order. Then we can solve the integral to find 



Im/ 



7T 



(d-3)* 



(d - 3) ( 2 d ~ 2 M ) d ' 3 E - [2M A - d Y^ E 2 + 



+aL 2 Pl E 3 ( (d - 3) ( 2 d -' z M 



ri-2*A 2 ( 2M i - d ) d ~ 3 E ) + 0(a 2 L Pl) 



The tunneling rate is therefore 

2vr 



r ~ exp 



(d - sy 



(d - 3) ( 2 d ~ 2 M ) d ~ 3 E - [2M 



E z + 



(4.6) 



+aL 2 Pl E 3 I (d - 3) ( 2 d -' z M 



^d-2n,A d-3 



2 ( 2M 4 " d ) d " 3 E ) + 0(a 2 L^; 



(4.7) 



For d = 4, the first and second expressions in the exponential exhibit non-thermal behavior 
that was discovered firstly in Ref. Q. So, the emission spectrum cannot be strictly thermal. 
In our situation, there is an additional term depending on the GUP parameter in first order 
and cannot be neglected once the black hole mass becomes comparable to the Planck mass. 
We now illustrate whether or not, the emission rates for the different modes of radiation 
during the evaporation are mutually related from a statistical viewpoint. Using ( [4.7[) , the 
emission rate for a first quantum with energy E\, gives 



lnF El 



2tt 



(d-3f 



(d-3)(2 d - 2 M) d ~ 3 E x 



2M 4 ~ d ) d ~ 3 E 2 + 



+aL Pl E-{ (d-3) 2 d - I M 



Similarly, the emission rate for a second quantum E2, takes the form 
2tt 



IiiTr, 



(d-3f 



(d-3) I 2 d ~ 2 (M - Ex) ) a ~ a E 2 - ( 2(M - E x ) 



1 

d-3 



4-d \ d " 3 z?2 



(4- 



+aL Pl E% ( (d - 3) ( 2 d - 2 (M - E x 



1 

d-3 



2 ( 2(M - Ex 



(4.9) 



- 15 - 



Correspondingly, the emission rate for a single quantum with the same total energy, E 
E\ + E2 , yields 



lnT 



2tt 



(d-3)* 



(d-3) 2 d ~ 2 Af M (£i+£ 2 ) 



1 

d-3 



2M 



4-d\ d -3 



(Ex + E 2 ) 2 + 



+aL I Pl (E l + E 2 f [(d-3) I 2 d ~ 2 M 



1 

d-3 



2 2M 



4— d \ d-3 



(4.10) 



It can be seen easily that these probabilities are correlated. On the other hand, the non-zero 
statistical correlation function is 



X (E 1 +E 2 ;E 1 ,E 2 



2itE 2 
(d-3f 



(d-3) 



2 d - 2 (M -Ex 



1 

d-3 



2 d ~ 2 M 



1 

d-3 



+ 



-En 



+aL Pl 4Ei ( 2M 



2M 4 



r4-d\ d -3 



2(M - E x 



+ Ei 2 a ~ z M 



d-2 71^4— d\ d -3 



+ 



1 

[2Ef + 2E\ + 3E x E 2 \-(d-3) (2^ 2 m) d ' 3 [3E\ + E\ + 3E X E 2 ] + 



+ (d-3)E% (2 d - 2 (M-E 1 ) 



2 d-2 M A-d\d-3 



2 d ~ 2 (M - E x ) 



1 - 




j d-3 


)] 







(4.11) 

Note that three first terms in square bracket in r.h.s. can be neglected for any reasonable 
values of quantities. Appearance of these terms originates from approximations used in 
the calculation of integral in equation (4.5). Therefore we find 



x(Ei + E 2 ; E\,E 2 



2iraL 2 P! E 2 



4£i 2M 



r4-d\ <*-3 



[2E 2 + 2E\ + 3£i E 2 ] 



(d-3) 2 

-(d - 3) (2 d ~ 2 M) ~ [3E 2 + El + 3ExE 2 ] + (d - 3)E% (2 d - 2 (M - Ex)) ~ + 

2 d - 2 (M-Exf 



2 d-2 M 4-d\l-3 



1 

d-3 



(4.12) 



This means that not only the probability of tunneling of two particles of energy E\ and E 2 
is not the same as probability of tunneling of one particle with their compound energies, 
Ex + E 2 , but also there are correlations between them. In fact, whenever one quantum of 
emission is radiated from the surface of the black hole horizon, the perturbations are created 
on the Planck scale that influence the second quantum of emission, and these perturbations 
cannot be neglected. Particularly, once the black hole mass becomes comparable with 
the Planck mass, one cannot neglect the effects of these perturbations. Therefore, as 
expected in Ref. [45], in this way the form of the amendments as back-reaction effects by 
incorporation of GUP effects, are sufficient by themselves to create correlation and therefore 
recovery of information. So, at least part of the information leaks out from the black hole 
as the non-thermal GUP correlations and this provides a basis to solve information loss 
problem. The question which arises here is time evolution of these correlations. This 
problem is currently under investigation. 
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5. Summary 

In this paper, generalization of the standard Hawking radiation via quantum tunneling 
through the smeared quantum horizon of the non-commutative Schwarzschild black hole 
has been studied. This study is based on the solution of the Eq. ( |3.5| ) in the context of 
coordinate coherent state noncommutativity and in a model universe with extra dimen- 
sions. The corrections of the emission rates for both massless and massive particles have 
been performed. We have shown that noncommutativity and extra-dimensional consider- 
ations are not enough to create correlations between emitted modes. That is, there are no 
correlations between the tunneling rates of the different modes in the black hole radiation 
spectrum at least at late-times and just to second order in E. However, spacetime non- 
commutativity is adequate by itself to preserve information due to the fact that in the NC 
framework, the black hole does not evaporate completely and this leads to the existence of 
a minimal non-zero mass (e.g., a Planck-sized remnant) in which black hole can reduce to 
it at final stage of evaporation. So, part of information might be preserved in this remnant. 
If one really believes the idea of stable black hole remnants due to the fact that there are 
some exact continuous global symmetries in the nature |]52|| , then, based on our findings 
(Sections || and |||), we should accept this possibility that a part of information stays in- 
side the black hole and can be retained by a stable Planck-sized remnant 9 . On the other 
hand, the generalized uncertainty principle of string theory or noncommutative quantum 
mechanics brings significant modifications to the Hawking radiation. As it is well-known, 
the GUP leads to a Planck-scale remnant at the endpoint of evaporation process. We have 
shown that the GUP as modification of the de Broglie relation in the quantum tunneling 
framework of the black hole evaporation, leads to correlation between emitted modes of 
evaporation. Hence in this setup, a part of the information leaks out of the black hole in 
the form of non-thermal GUP correlations, and the other part reside in the remnant. These 
features have the potential to answer some questions regarding the black hole information 
loss paradox and provide a more realistic background for treating the black hole evapo- 
ration process, especially in its final stages of evaporation. We stress that time-evolution 
of GUP-induced correlations is an important issue which needs further investigation and 
probably shed more light on information loss problem. 
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The imaginary part of the action (Eq. ( |3.2C| )) with some arbitrary number of dimensions 
9 Of course, this assumption is acceptable if information is really conserved in our world. 
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can be calculated as follows 



d = 4 — ► Iml = 2^G 4 e~ — El -4E—^ — 4— +2E— [ — ) e~ 
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